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^^ ■ Abstract. The Generator Coordinate Method (GCM) in the Gaussian Overlap 

Approximation (GOA) is applied to a description of the nuclear quadrupole collective 
states. The full five-dimensional quadrupole tensor is used as a set of the generator 
coordinates. The integral Hill- Wheeler equation is reduced to a differential equation 
by using the Fourier transforms of the overlap and energy kernels. The differential 
y ' Bohr Hamiltonian obtained this way is compared with that derived by the usual 

(-H . approach to the collective Hamiltonian in the GOA which does contain an additional 

approximation. The method of calculating the quantities which determine the Bohr 
Hamiltonian from the set of deformation-dependent intrinsic states is demonstrated. 
^ , In particular, it appears that the moments of inertia at the quadrupole rotations are 

Cn ' of the type of that of Yoccoz. 
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1. Introduction 

The Bohr Hamiltonian (cf a recent review [Ij) still remains one of the most effective 
tools to describe the lowest collective quadrupole excitations of the even-even nuclei. 
This is in spite of a long history of the collective model dated from the fifties of the 
twentieth century j2]. The rudiments of the model can be met even earlier [3]. Nowadays 
the Bohr Hamiltonian in a general form proposed in [H |5] is considered. It is based on 
the correspondence principle, that is, on the assumption that the collective quadrupole 
motion has its classical counterpart. The quantum Bohr Hamiltonian is constructed 
from the most general classical Hamiltonian quadratic in velocities by the Podolsky- 
Pauli quantization prescription [HI [TJ [8]. The quantum kinetic energy operator is then 
the Laplace-Beltrami operator in the Riemannian space with the metric tensor equal to 
the classical inertial matrix. Derivations of the Bohr Hamiltonian from the microscopic 
many-body theories of nuclei used so far, which are based on methods of the type of the 
Time-Dependent Hartree-Fock method, do pass through their classical stage too. The 
collective motion is then represented by a wave packet, which keeps the classical law of 
motion, requantized afterwards (cf e.g. [9l[l]). 

A possibility of a purely quantum description of collective excitations is created 
by the Generator Coordinate Method (GCM) ^ [11] (see also [I2|). The method 
leads originally to an integral eigenvalue equation (the Hill- Wheeler equation). The 
integral equation can be approximated by a differential equation of the Schodinger type 
[13j. It is especially easy to do when the Gaussian Overlap Approximation (GOA) is 
made for overlap and energy kernels in the Hill- Wheeler equation [TT| [T^ [T5] . The 
method is usually demonstrated for a single generator coordinate. But the quadrupole 
motion has five degrees of freedom and thus five generator coordinates are required. 
The mult i- dimensional GCM was already considered by several authors |16[ [T7t [T8| [T9] 
(see also the review article [20]). Obviously, the GCM, being one of the contemporary 
tools for investigations of various collective phenomena and for the restoration of broken 
symmetries, has been applied to various problems in the nuclear physics apart from the 
quadrupole collective motion. The two-dimensional Hill- Wheeler equation has been used 
to investigate the quadrupole /3 and 7 vibrations [2 H [22 | |23 ] . Before that the differential 
equations coming from the GOA have been applied to these vibrations |2^ [25l [26l |27] . 
The GOA has been used to estimate the ground-state correlations associated with the 
quadrupole collective modes [28]. But the rotational correction to the ground-state 
energy has been obtained through the projection onto the good-angular-momentum 
state. The five- dimensional collective quadrupole Hamiltonian based on the GOA has 
been generated from the microscopic theory with the Gogny effective forces and used in 
comprehensive calculations of the 2+ nuclear excitations [29j. However, the moments of 
inertia obtained by the Thouless-Valatin prescription [30] instead of those derived from 
the GOA are used in the Bohr Hamiltonian. In spite of these developments, it seems 
there is no consistent derivation of the full five- dimensional Bohr Hamiltonian using the 
GCM with the GOA yet. 
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The present study is devoted to application of the multi-dimensional GCM in the 
GOA to the five- dimensional quadrupole collective motion of the even-even nuclei. The 
quadrupole deformation tensor plays a role of a set of the five generator coordinates. 
We aim at deriving the Bohr Hamiltonian in a purely quantum way. The Hamiltonian 
obtained in such a way is a bit more general than that based on the correspondence 
principle (cf [Ij). In order to transform the Hill- Wheeler integral eigenvalue equation 
into the Schrodinger differential equation with the Bohr Hamiltonian we introduce, after 
Onishi and Une [IT], new quadrupole variables which diagonalise the exponent of the 
Gaussian overlap. Using these new variables we perform the Fourier analysis of the 
energy kernel. The kernel turns out to be a local distribution containing the Dirac delta 
function and its derivatives of the second order. It allows us to reduce the integral 
equation to a differential one. The Hamiltonian of the rigid rotor was derived earlier 
using the same method [31]. We compare the method of the Fourier analysis with the 
method used usually in deriving the collective Hamiltonian from the GCM [12| [201 [T9] . 
The latter method contains an additional approximation apart from the GOA. 

We formulate the Hill- Wheeler theory for the quadrupole collective excitations 
in section [2l Section |3] is devoted to the Gaussian Overlap Approximation for the 
quadrupole generator coordinates. The differential equation equivalent to the Hill- 
Wheeler equation in the GOA is derived and compared with the usual collective approach 
to the GCM in section 14.11 Formulae for calculating the quantities determining the 
overlaps within the GOA and the Bohr Hamiltonian from the deformation-dependent 
intrinsic state are given in section [5l We draw conclusions from the present study in 
section [6l In Appendix A the semi- Cartesian Wigner functions are introduced, which 



are convenient to use in the present study. The structure of matrices being isotropic 



functions of the quadrupole variables is discussed in [Appendix B 



2. The Hill- Wheeler theory of the quadrupole collective excitations 

2.1. The quadrupole generating states 

Let the considered nucleus with Z protons and A^ neutrons be described microscopically 
by the many -body Hamiltonian H which is invariant under the 0'^{3) group 
of transformations i.e. the superposition of time reversal T and the orthogonal 
transformations (the rotations R{u) and the space inversion P) in the physical three- 
dimensional space. In order to describe its quadrupole collective excitations by means of 
the GCM we start with states parametrised by the two deformation parameters do and 
d2 which describe the quadrupole deformation with respect to a system of principal axes. 
For instance, the intrinsic components of the mass quadrupole moment in the deformed 
ground state can serve as such parameters. Parameters c/q and ^2 define the deformation 
(e.g., a shape) up to the group of five permutations of the principal axes (cf [T]). One of 
the principal axes is the symmetry axis when one of the three following conditions for 
the deformation parameters is fulfilled: d2 = 0, d2 = ±-\/3rfo- The conditions represent 
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the three straight hnes (or the six rays) intersecting each other at point (io = 0, (^2 = 
on plane (^0^2) of the deformation parameters. The rays divide the plane into the six 
sectors. The values of d^ and 62 located in one of the sectors describe the all possible 
deformations (shapes) at the established names of principal axes. The standard choice 
of the sector is the following: 

< 4 < VSrfo < +00, (2.1) 

which means that the symmetry axes are z and y for (i2 = and ^2 = Vorfo? respectively. 
Here, we assume that the ranges of values of do and d2 are given also by (12. ip . We 
denote for short the set (rfo, ^2) as d. Let |(/>((i)) be an approximated ground state 
of the considered nucleus with a given deformation d. State |0((i)) is supposed to be 
the ground state of the mean-field Hamiltonian, either the phenomenological or the 
self-consistent one with constraints for the quadrupole moment. For all possible values 
of d states |0((i)) form a set of states called the intrinsic states in the following. The 
normalization of the intrinsic state is assumed to be: 

{m\4^{d)) = l (2.2) 

for every d. 

In general, the state |0(c?)) is not an irreducible representation of 0"^(3) because 
the rotational symmetry is usually broken spontaneously in the mean field. However, it 
is assumed that \(t){d)) is still invariant under the point D^^ group j32] with the intrinsic 
axes (x, y, z) as its symmetry axes and we have: 

f\(t>{d)) = \m), (2.3a) 

e-^--^"|0(c/)) = |0(c^)), {2.3b) 

P\cP{d)) = |0(rf)), (2.3c) 

where J^ for u = x,y,z are the intrinsic components of the total angular momentum 
operator. Equation (I2.3ap can be fulfilled only for even Z and N and, as a matter of fact, 
it corresponds to the assumption that we will consider only the even-even nuclei here. 
Although the Djf^ symmetry looks natural, the parity invariance (12.3d) will not be used 
explicitly in further considerations. It is because the quadrupole deformation tensor 
itself is of the positive parity. Since the differentiation with respect to the deformation 
parameters d commutes with the Dj^ generators the states 

\Md)) = ^m)), I'^'^K^)) = ^|^I0(^)), ••• (2.4) 



for k, 1 = 0, 2 have also the D^^ symmetry and fulfil (12.3fflp . ( 12.3&I) and ( 12.3cp . 



It follows from time reversal invariance that the matrix elements obey relations 

(0(rf)|0(rfO) = {<l>{dM{d')y = {<i>{d'M{d)), 

{<P{d)\H\<P{d')) = {<P{d)\H\<P{d')y = {<P{d')\H\<P{d)), (2.5) 

i.e. they are real and symmetric with respect to d and d'. On the other hand, 
matrix elements {(f){d)\Ju\(/){d')) and (0((i)| J„^|0(rf')) are imaginary and antisymmetric, 
because the angular momentum operators Ju are time odd. 
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It follows from the D2 symmetry (12.3&P of |0(c?)) and |0i((i)) that, in particular, 
the matrix elements of type {(j)i{d)\Ju\(t>{d)) , {(j)i{d)\JuH\(f){d)), and {(f){d)\JuJu'\(p{d)), 
{(f){d)\JuJu'H\4>{d)) for u ^ u' all vanish. 

Since \(f){d)) does not possess the rotational symmetry, the state rotated by the three 
Euler angles u = {ui,U2, W3) from the ranges < Wi < 27r, < a;2 < tt, < Wa < 27r 

\^{d,uj)) = R{uj)\(f){d)) (2.6) 

is an equally good ground state of the mean-field Hamiltonian. The unitarity of R{u) 
assures us that it is also normalised to 1. Operator R{u)) is time-even and thus the 
matrix elements between states \^{d,uj)) and \^{d',u')) corresponding to these of ( 12. 5 P 
are also real and symmetric with respect to sets of variables {d,u) and {d',u'). States 
(12. 6p will play the role of the generating states and the five variables do, d2, Ui, lo^ and 
bjj, can be used as the generator coordinates for the quadrupole collective excitations. 
As a matter of fact, other generator coordinates, namely the laboratory quadrupole 
deformations defined below will be used. 

2.2. Quadrupole coordinates 

We introduce a new set of the five variables, namely a^, // = —2, . . . , +2, by the relation 

a^ = doD^o{uj) + d2Df,2iuj), (2.7) 

where the semi- Cartesian Wigner functions D^o and -D^2, defined by ( lA.ll) . are 
combinations of the usual Wigner functions T>'tj^. They play the role of the laboratory 
deformation parameters. The Jacobian of transformation (12. 7p is equal to (cf [Ij) 

W{d,uj) = \d2{Ml-dl)\ sin UJ2 (2.8) 

and, therefore, the transformation is reversible in the entire assumed open ranges of 
variables d and u. It is seen from (12. 7p that a^ are the components of an electric (i.e. 
with positive parity), real (i.e. a* = (— l)'^ci;_^) quadrupole tensor ex in the three- 
dimensional physical space. The complex conjugate components are denoted with the 
superscript, e.g. 

«^ = «; = (-i)'^«-M- (2.9) 

The frame of axes with orientation w = is the frame of principal axes (intrinsic 
frame) of a, because the intrinsic components of ck are equal to ao{u) = 0) = do, 
a2{uj = 0) = a_2(co' = 0) = d2/V^ and ai{uj = 0) = a_i(w = 0) = 0. 

Instead of the complex spherical components of the tensor, the five truly real 
coordinates a^, k = 0, 2, x, y, z, namely 

Ofc = Dl{uj = 0)«^ = Cj^a^ (2.10) 

can be used, where the Einstein contraction rule is applied to summing of the Greek 
upper and lower indeces fi = —2, . . . , +2 of the spherical components. The rule is 
not applied to the Latin Cartesian indeces. These coordinates can be treated as the 
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Cartesian components of vector a in the five- dimensional space. This is because the 
scalar product of the two quadrupole tensors, ex and (3 is 

ocj3 = J](-l)^a^/3_^ = a^/3^ = J^ a^hk = a ■ b. (2.11) 

M A; 

In connexion with (12. lip the volume element in the space of the quadrupole coordinates 
is: 

dn(a) = Ilfcdafc = 14(3^0 — dl)\ sinc<;2d(iod(i2dco'idco'2dw3. (2.12) 

All functions of deformation d and the Euler angles, u, can be treated as functions 
of either tensor a or vector a. Here, we prefer to use the complex variables, a^, for 
convenience. However, we have in mind that the integration variables are always real: 
either a^ or d and u as is seen in (I2.12p . The relation between the derivatives with 
respect to a^ and the derivatives with respect to the Euler angles and the intrinsic 
components do, 4 reads |4H [T]: 



daa ^ dttk 



odn 



->dd. 


+ 'D^H 


' L 

24 


"'-'^'<"'2i. 


Ly{u] 



- iD^,{uj)—L,{uj) - iD>^{uj)—Ly{uj), (2.13) 

where 

4 = -2(v^4 + 4), dy = -{^/?,dQ- d2), 4 = 4 (2.14) 

and Lx{u), Ly{u)), Lz{uj) are differential operators with respect to the Euler angles 
and are interpreted as the intrinsic components of the drift angular momentum of the 
rotation of the intrinsic frame (eq. (2.15) in ref. pOj). 

2. 3. Trial states for the collective quadrupole excitations 

We denote the generating states (12. 6 P as |$(q:)). Then trial states for the quadrupole 
collective motion are taken in the form: 

\^iJ[ip]) = I ip{ct)\^{a))dn{<x). (2.15) 

The five components of the quadrupole tensor a play the role of the generator 
coordinates and function <f{oc) is the weight function. The variational principle leads to 
the Hill- Wheeler integral equation of the form 

f[n{a, a') - EI{a, (x')]^{(x')dQ{(x') = (2.16) 

for function ip of the generator coordinates. The equation is determined by the two real 
symmetric kernels: the overlap kernel 

X(a, a') = (*(")!*(«')) (2-17) 
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and the energy kernel 

H{a,cx') = {^{cx)\H\^{cx')). (2.18) 

The overlap kernel fl2.17p is non-diagonal and equation ( I2.16P for the weight functions 
constitutes a non-orthogonal eigenvalue equation. 

3. The Gaussian Overlap Approximation 

3.1. Gaussian approximation for the kernels 

We expect that the overlap kernel drops quickly from one to zero when differences 
In = (y■^l — ttu of the components of tensors ot and a' increase from zero. Therefore, 
the logarithm of the kernel rather than the kernel itself as a function of 7 can be 
approximated by the power series: 

ini(/3 + It, /3 - \i) = -Ig'^iOh.^. + • • • , (3.1) 

where f3 = |(ck + Q;') and the Einstein contraction rule is applied to upper and lower 
indeces fi, v. When in expansion (13. ip we keep only the term quadratic in 7 for an 
arbitrary /3, we have the local Gaussian approximation for the overlap kernel: 

X(/3 + |7,/3 - h) = exp {—g^'^ia)!,!,). (3.2) 

Matrix g is real, symmetric and positive definite. We use a similar approximation for 
quotient 

Matrix h is real and symmetric. 

Both matrices, (7^^ and /i^j^, are the symmetric quadrupole bitensors in the 



interpretation of Appendix B whereas v a scalar. We shall refer to matrix g as the metric 
tensor. Throughout the paper the upper and lower indices of matrices are connected 
with the complex conjugation of the spherical tensors according to (12. 9p . and not with 
the contravariant and covariant components of tensors in the Riemannian space. A 
separate symbol is used for the inverse metric tensor, namely g~^ = f. 

3.2. Eigenvalue equation for the Gaussian overlap kernel 

To find the eigenfunctions of the Gaussian overlap kernel ( 13. 2p it is convenient to reduce 
the exponent of the Gaussian function to the sum of squares of five variables. It is a 
simple thing to do for the single generator coordinate |33] . It is not so simple in the 
case of a set of a few variables. Here, we adopt the defininion of the new variables given 
by Onishi and Une [I^. Namely, we introduce the five new real variables tk {k =0, 2, 
X, y, z) through a line integral in the five- dimensional space of vectors a of (I2.10p 

tk{a)= I ^rfcKa')da;. (3.4) 

"^ I 



GOA for the quadrupole states 8 

The lower integration limit is omitted in the notation of the integral because it is not 
relevant here. Matrix r is taken to be the positive-definite symmetric square root matrix 
of g defined by ( IB. 81) with m = g and m^/^ = r. It fulfils relation ( IB. 71) . namely 

gki{a) = ^rfci(a)rj,(a), (3.5) 

i 

where Qki are the Cartesian components of g defined according to (lB.4p . The variables 
tfc are defined well if the line integral does not depend on the path of integration. It is 
so for the irrotational field rki{a) i.e. for 

drki _ drki_ ^ ^ ,^^. 

dai dai 
which we assume for use of further considerations. However, we do not know in advance 
whether it is really fulfilled. 

It is convenient to operate with the complex linear combinations r^ of tk defined 
according to ( I2.10p . Then we have 

dr^(a) = ^ C^kC^rki{a)dai = r'j[{a)dau, (3.7) 

k,l 

where r'^ = r^ = r'^ ^ are the corresponding spherical components of r (the sequence 
of the superscript and subscript does not need to be fixed for a symmetric matrix). 
Relations (13. 5p and (13. 6p have the following counterparts for the spherical components: 

) -Y,rH{a)ra{a)\\ =0, (3.8) 

= 0. (3.9) 



g^^{a) - r^{cx)rxu{a) = 2^ 


C^kCul Qk 


k,i 


\ \ 


k,Li 


(drki drki 
\ dai dai 



The Jacobian of transformation r = t{oi.) is equal to ^Jg{ol) > 0, where g = det (g). 
The transformation is then reversible. The reversed transformation is denoted as 
ot = cx.{t). The volume element in the space of r is 



dfi(T) = ^/g{a)dn{a). (3.10) 

When we assume that the values of 7^ = a^ — a' are small, what is consistent with 
the Gaussian approximation, we have from ( 13. 7p : 

^, = t,-t'^ = r; (l(a + a')) (a. - a',) = r;(/3)7., (3.11a) 

0,=ki^, + r'^) = r,{(3), (3.116) 

where r^ = t^{cx'). 

Using the new variables, r, it is easy to solve the eigenvalue equation 

y"x(«(T),a(T'))x(6,r')dfi(r') = j(e)x(e,r). (3.12) 

for the Gaussian overlap. The wave tensor e stands for a set of the five quantum 
numbers, e^ {fi = —2, . . . , 2), which label eigenf unctions x(e, r) and eigenvalues j(e). 
The normalised eigenfunction reads: 

x(e,T) = — -^exp(ie-r). (3.13) 
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It fulfils the following orthogonality and completeness conditions: 

j x*(e, T)x(e', T)dfi(r) = 5,{e - e'), (3.14a) 

j x*(e, T)x(e, T')dfi(e) = 5,{t - r'). (3.146) 

The tensor properties of quantum numbers e^ are clearly seen from (I3.13p . The 
eigenvalue is equal to: 

j(6) = (27r)5/2exp(-i6-e). (3.15) 

3. 3. The orthogonal Hill- Wheeler equation in GOA 

We know that the Hill-Wheeler equation (I2.16P can be formally reduced to the 
orthogonal eigenvalue equation (cf [12j) for wave-function il){(y.) defined as follows: 

^{a)= n{a,a')y^{a')dQ{a'), (3.16) 

where 

7^(«, = 1 Xie, T{a))j'/\e)x*ie, T(^))dfi(e) (3.17) 

is the square root kernel. The name of TZ is justified by the relation 

f n{a, 07^(^, a')dfi(r(0) = X(a, a'). (3.18) 

Equation fl3.16p can be formally reversed and the weight function, if can be expressed 
in the form: 

^(«) = JATia, 0^(Odfi(T(0), (3.19) 

where the narrowing kernel jTT] 

A/'(«,0 = y"x(e,T(a))j-i/2(e)x*(e,T(0)d^(6) (3.20) 

is introduced, which fulfils the following relation: 

J 7^(a, OAr(|, «')dfi(T(0) = 55(T(a) - T(a')). (3.21) 

We define the narrowed overlap kernel 
Jita = y"Ar*(a,0X(a(T),a'(T))Ar(«',Odf^(T(a))dfi(T(a')) 

= 55(t(0-t(O), (3.22) 

which takes the diagonal form. 

To obtain the integral eigenvalue equation for function ip we substitute (I3.19P for 
the weight function, ip in the Hill- Wheeler equation (I2.16p . multiply it by Af{ct,^) and 
integrate with respect to T{a.). Then we get 

J /C(^, O^(Odn(T(0)) = E^{^), (3.23) 
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where 

X:(^, ^') = I U*{a, ^)nicx, a')^f{a', Odf^(T(a))df](r(a')) (3-24) 

is the narrowed energy kernel. 

3.4- The Fourier analysis of the Gaussian energy kernel 

The problem with equation ( 13.231) . and kernels ( 13.201) and ( I3.24p is that eigenvalues j(e) 
of the Gaussian overlap kernel tend to zero when e ■ e tends to infinity. It means that 
the integral in (I3.20p is divergent and we have the "ultra-violet catastrophe" for A/". 
To make the integral convergent we can introduce a cut-off for e • e. Instead, we can 
use operational methods of the distribution theory and treat integrals being the Fourier 
transforms of the non-integrable functions as convergent due to oscillations with the 
infinite frequencies of the integrand and express them by the Dirac delta function and 
its derivatives, namely 

1 r+oo 

— j f{k)e''^dk = f{-id/dx)6{x). (3.25) 

Equation (I3.25P follows from the formula for the derivatives of the delta function 
represented in the form of the Fourier integral, namely: 

-i— 5(x) = — / fcV^^dfc (3.26) 

\ dxj 27r J_^ 

for n = 0, 1, .... 

In order to express the Gaussian energy kernel in terms of the delta function and 
its derivatives we introduce the Fourier transform, J-", of "H in the following way: 

^(e, e') = j x*{e, t{(x))H{cx, cx')x{e\ T(a'))d(](T(a))df](T(a')). (3.27) 

Substituting (Q, (l33T^ . (l3TT&|) and (l37[3|) into (131271) we have 



X exp (-ia ■ e{f3)) exp (-i^(/3, 7) ■ ^(/3, 7) - i^ ■ ^(/3, 7))dfi(^(/3))dn(^(/3, 7)), (3-28) 

where 

a.x{l3) = h''''i(3)n^.{(3)nxui(3) (3.29) 

and e = |(e + e'), cr = e — e', and n^j, = iXl ^/g){d^/dr^^) is the matrix reversed to 
r, i.e. Uyxr^^ = 6^. The Gaussian integrals with respect to t? can be easily calculated 
and we obtain 

•^(e, e') = ^ J[v{f3) + lie^ex - (5^)a^(/3)] exp {-ia ■ 0{f3))dn{0{f3)). (3.30) 

The narrowed energy kernel (13.241) expressed in terms of the Fourier transform 
(13:271) reads: 

m, ^') = j 3-''\e)3-''\e')x{e. ^)J^ie, e')x*{e\ r)d(](e)d(](e'). (3.31) 
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Inserting (IXT^ . (EED and ^M) into flg:^]) we have 

m, O = J^ J exp (|o- ■ ^) exp {icT . {e{rj) - e{f3))dn{cT) 

X [v{f3) + lie^e, - 5l)alm exp {is ■ ^(77, C))dfi(£)dfi(0(^)), (3.32) 

where ^7 = i(^ + ^') and C = ^ ~ ^'- The integrals with respect to e and cr in fl3.32p are 
the fivefold integrals of type (I3.25P and thus can be expressed in terms of the partial 
derivatives of the five-dimensional Dirac 5 function in the following way: 

^ ^[v{(3) + \{e-e^ - ^Dalm exp (is ■ ^(77, C))df^(s) 



(27r)5 



d 



d 



SI <(/3) 



^ j exp (|o- ■ a) exp (io- ■ (^(ry) - 0(/3))dfi(^ 



^5(^(^7,0), 



(3.33) 



(2vr) 
= exp 



— 1- 



55(0(r7)-0(/3)). 



(3.34) 



After inserting (13.331) and (I3.34p into (I3.32p the integration with respect to 0{/3) is easily 
performed and finally the narrowed energy kernel takes the form: 

/C(eO = exp(-|A(0(r7))) 



X 



v{r)) 



9^«(r7,C)9^^(r7,C). 
where f = g^^ = n ■ n is the matrix reversed to g, and 

d d 



H^ivX)), 



(3.35) 



Mr) 



Oth ^r^^ 



(3.36) 



The kernel is thus the distribution dependent on the delta function of 'J9(^,^') = 
r(^) — t(^') and its partial derivatives of the second order. We rewrite it in a shorter 
form: 

d d 



/c(^,r 

where 



vil) - ia^(r7) 



d^^{riX)d^^ir]X) 



S,{^{riX)), 



(3.37) 



a/iuirj) = exp (-|A(0(r7)))a^^(r7), (3.38a) 

v'{r)) = exp (-iA(0(r7))) [v{ri) - lh^''{ri)UM] = exp (-|A(0(r7)))t;'(r7). (3.386) 

4. Differential eigenvalue equations 

4.I. The differential Hill- Wheeler equation 



The integral Hill- Wheeler equation (I3.23P is easily converted into a differential equation 
using the kernel in the form of (I3.35p . Indeed, inserting (I3.35P into (I3.23P we have: 



^(^') 



dT-{^') dr^i^') 



55(T(0-r(O)dn(r(O) 



Ei;{^). 



(4.1) 
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Integrating by parts we obtain from (14.1 p the following differential eigenvalue equation: 

m) = Em). (4.2) 



Id d - ' 



where 



Coming back to the original variables ^ in (14. 2 p we transform the derivatives with 
respect to 'r^(^) into the derivatives with respect to ^^ and we obtain the eigenvalue 
equation 

H{i)m) = E^{^) (4.4) 

for the Bohr Hamiltonian 

with the inverse inertial bitensor 

V(6 = n,.{^)n,,{^)-a''\^) (4.6) 

and the potential 

^(e=exp(-|A(5))„'K)---^-(v5(«n„„(e"..(fl^) . (4.7) 

The Laplacian (I3.36P takes the form: 

The relation 

J- (v/M^..(0) = 0, (4.9) 

which comes from (13. 6p . has been used to obtain (14. Sp . (14. 7p and (14. 8p . 

Equation (14. 4p is equivalent to the Hill- Wheeler equation (I3.23P with the Gaussian 
plus quadratic approximation (13. 2p and (13. 3p for the kernels. No other approximation is 
made. However, one can say that the equivalence is formal only because the quantities 
(14. 6 P and (14. 7p can be determined merely in an approximation when expanding the 
exponential operator. For instance, we expand (13.38 ap : 

a,.(0 = «Mv(0 - l^^WM) + \ (^) A(OA(Oa,.(|) + • • • (4.10) 

However, we have no guarantee whether such an approximation procedure is convergent. 
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4-2. The collective Hamiltonian 

Another method of deriving the Bohr Hamiltonian, which passes by the integral 
equation (I3.23P and starts directly from the variational principle, is more often used 
(cf [T2| [T9| |20]). We refer to it as the collective approach to the GCM. It consists 
in representing the energy mean value in the form of the expectation value of a local 
differential operator, say if, Hermitian with weight w in the space of the generator 
coordinates, namely: 



i)*{oi)H{oi)i){oi)w{cx)dn{oi). (4.11) 



In the GOA the square root kernel fl3.17p is equal to 

^(«, i) = i^-) exp (- (T(a) - t(^)) • {Tic,) - t(0)) (4.12) 

and then, according to (I3.18P and (13. 3 P the energy kernel can be presented in the 
following form 

nia,a') = Jnicx,^M(3) - i/^'^^(/3)7,.7.]7^(^,a')^^(OdO(0 

= j n{cx,i)[v{l3) - ^a'^'^(/3)^^(a,a')^.(a,a')]7^(^,a')dfi(r(0). (4.13) 

An additional approximation should be made to achieve the form (14. lip for the 
energy functional. Namely, we assume that functions f (/3) and a^^{(3) do not change 
much within the range of /3 = |(ck + ex') in the vicinity of ^ where both Gaussian 
functions, TZ{cn, ^) and 7?.(^, ex') take still big values. In other words we assume that 
v{f3) and a^'^(/3) are slowly varying functions in the vicinity of ^, and we can expand 
them in the power series around $,: 

vm ^ viO + ^(/3, - U + ||^(/3m - e.)(/3. - e.) + . . . (4.14a) 

a^''{(3) ^ a'^'^i^) + ^^^(/3. - Q + i^^^(/3, - e.)(/3A - 6) + • • ■ (4.14&) 



9e« ^"^ "^^ 2 Q^^g^^ 



A 



in (I4.13p . Here, we will exploit the lowest approximation and take the zero-order terms 
only. Since the Gaussian functions with monomial factors in front are equal to the 
derivatives of the Gaussian functions themselves, the (x- and a'-dependence of the 
integrand in (I4.13P appear only in functions TZ{cx., ^) and 7^(^, a') and their derivatives 
with respect to ,^^. Then, the integration with respect a and tx' in (14. lip can be 
performed giving the collective wave-functions ip*{$,) and ip{$,) of (I3.16p . Finally, we can 
obtain the energy mean value in the form (14. lip performing integrations by parts. The 
collective Hamiltonian in the space of coordinates t(^) reads 
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with a^y given by fl3.29p and 



8ar^(0ar.(0 



= .-(0-1 ^'"^-^^^ . (4.16) 

Coming back to variables $, we end up with the collective Bohr Haniiltonian 

Hermitian with weight w{$,) = \/ g{$,), where 

A,M)=f,M)U{iW\i) (4.18) 



v{i) = vio - \h^^{i)f,M 



^ ^ ^ ( ^)n,M)n..{if-^^ 1 . (4.19) 



and 



We see that the Bohr Haniiltonians H and H have the same form which is identical 
with that discussed in [T] for the quantum collective models. The only difference between 
them is that they are determined by different inverse inertial bitensors and different 
collective potentials, A^i, and A^^, and V and V, respectively. For exp(— |A) k, 1 
the quantities with bar transform into the corresponding quantities without bar. Such 
an approximation is good when a^y{$,) and v{$,) are slowly varying functions of their 
argument. It is in accordance with the assumption made when deriving Hamiltonian H 

of dnn). 

4-3. The zero-point energy 

In the original approach to the Bohr Hamiltonian, there is an ambiguity related to the 
correspondence principle. The Podolsky-Pauli quantization prescription tells us how to 
quantize the kinetic energy; however, in general, the kinetic energy operator is given up 
to an additive arbitrary function of the deformation [8J. In consequence, the potential 
energy which enters the Bohr Hamiltonian operator needs not to be identical with the 
classical one, as has been often assumed in the past (cf [Ij). In the present study it 
would correspond to the collective potential V{d) = v{d) = {(f){d)\H\(f){d)) equal to 
the static ground-state energy term . On the other hand, it is obvious that dynamical 
correlations in the ground state of quantum systems should appear. Various effects 
of the ground-state correlations have been investigated for a long time (see, e.g., a 
review [M])- In the theory of the collective excitations, the difference between the static 
energy and the collective potential is interpreted as the zero-point energy associated 
with a given collective mode. The GCM gives a definite result for the zero-point energy 
correction to the static ground-state energy. The zero-point energies associated with the 
quadrupole vibrational modes have been already estimated in the frame of the GOA, 
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e.g., |36[ Eni Ell 1221 EZ]- Those associated with the rotational modes were investigated 
rather by means of the angular- moment urn projection technique [22| |37] . 

In the present five-dimensional GCM with the set of generator coordinates forming 
the quadrupole tensor, all the quadrupole modes, both vibrational and rotational are 
treated on an equal footing. The collective potential of fl4.19p contains the zero-point 
energy correction, f(^) — V^(^), associated with all the five quadrupole modes. We 
see from f l4.19p that the correction is composed of two terms. The first, |f ■ h, is 
connected with the quadratic term in the energy kernel expansion ( 13. 3 p and manifests 
the dynamical effects. The second one appears only when the coefficients in the Gaussian 
expansions (13. 2p and ( 13. 3p . g^u{(3) and h^^{f3), depend essentially on f3 (the case of 
the local Gaussian approximation). The additional zero-point energy corrections come 
from expanding the exponential Laplace operator in ( 13.38 ap and (13.386P as discussed in 



section 14.11 Then, additional terms appear in both, the inverse inertial bitensor, A^,^, 
and the potential, V, in the Bohr Hamiltonian H of (14. 5p . It seems that corrections 
of that type to the zero-point energy have been taken into account in [21] for the one- 
dimensional Bohr Hamiltonian. In our case of the five-dimensional quadrupole modes a 
counterpart of the collective potential with similar corrections would have perhaps the 
form: 

m) = (1 - |A(«)),/(«) - --^^ (^^MTk^MM^)^) . (4.20) 

The standard procedure of deriving the collective Hamiltonian, shown in [T2 | [T9 | |20] . 
consists in expanding quotient 'H{a,a')/X{a.,(x') of ( 13. 3p (the expression inside the 
square brackets in (I4.13P ) as a function of ck and ex.' in the power series around an 
arbitrary point $, up to the second order. The actual meaning of the approximation is 
not very transparent. Certainly, it constitutes a stronger approximation than the usual 
quadratic approximation of (13. 3p in the GOA. It seems that a counterpart of such an 
approximation in the present approach is expanding the static potential ( I4.14ap up to 
the second order. Then, the resulting collective potential would have the form: 

V{^) = (1 - |A(e))^(0 - |/^^'^(0/m.(0 
1 1 d 



Sv^M^^p 



4.4- Matrix elements of observables 



v/^n,.(On..(0^^) . (4.21) 



Solving the Bohr equation (14. 4p (or that with Hamiltonian (I4.17P ) provides us with the 
energies, Ej^i, and the collective wave-functions, ipj^i, of the collective states 

|vl/^) = y"A/'(«,OV^^(OI'^(«))V^dfi(^)\^Rdfi(«), (4.22) 

where A^ stands for a set of quantum numbers of the state. The states form an 
orthonormal set, i.e., (\l/7v|^I/Ar') = Snn' when 

J rN{^)^N'{^)V9ii)dni^) = 5nn'. (4.23) 
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To know other characteristics of the collective states we should calculate the matrix 
elements of the other observables. In the GOA we approximate the matrix element of 
observable Q within states |$(q:)) in the way similar to (13.31) : 

($(a)|g|$(«')) = g(/3) + ign/3)7. - \q'''mi,i. + ... (4.24) 

Proceeding further in the way similar to that of sections 13.41 and 14.11 for H we obtain 



2^^^^v.v..-.v.. v.v^.-.v..^^^ 



9{^)Qi^ ^7V'(0df^(0, (4.25) 

where Q/^i/ and Q are given by (14. 6 P and (14. 7p when h"'^ and v are replaced with g"^ 
and q, respectively, whereas 

g,(0 = n,^{^) exp (-|A(0) r (0<(^)] , (4.26) 

and ^jv are eigenf unctions of Hamiltonian H. 

When q, n'^qu and n'^q^xn'^ are slowly varying functions of $, the exponential 
Laplacian operator can be replaced by the unity and we have the expression similar 
in form with (I4.25p (cf |25]). namely: 



(vl>;v|Q|^7V') = hlri^) 



^|^v/^gM.(0|;; 



'^^ W)QM) + vW)Q,i^)-^ 



2 ^^^^ V .V..-..V.. v.v^.^.v..^^^ 

+ ^/MQiO] i^N'im^m, (4.27) 

where Qfj.u, Qfi and Q are given by the corresponding formulae for Qfj,u, Qfi and Q with 
the exponential operator put equal to unity, and tpN are this time eigenfunctions of 
Hamiltonian H. 

5. The deformation dependence of the quadrupole collective Hamiltonian 

The original object in the present approach to the description of the quadrupole 
collective states in even-even nuclei is a set of the many-body intrinsic states \4>{d)) 
parametrised by the two quadrupole deformations d = {do,d2). Knowing \4>{d)) we can 
express the overlap and energy kernels, the two basic quantities in the Hill- Wheeler 
equation, in terms of overlaps and matrix elements of H and J^ within states \4>{d)) 
themselves and their derivatives (12.40 . All these overlaps and matrix elements are 
functions of the deformation. Having the kernels as functions of the deformation 
we calculate the collective potential and inertial functions which determine the Bohr 
Hamiltonian. 
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5.1. The overlap and energy kernels as functions of deformation 

The overlap kernel X(a, a.') is real, symmetric and normalised to unity for the all values 
of variables a^, namely 

X(a, a) = ($(a)|<l>(a)) = 1. (5.1a) 

Differentiating twice (15. lap we have the following conditions: 

($(a)|<l>^(a)) = 0, {5.1b) 

(<l>(a)|<l>^'^(a)) + ($^ (a) !$''(«)) = 0, (5.1c) 

where 

\^^{c.)) = ^|«1>(«)), |<|.^-(a)) = ^^l$H) (5.2) 

are the states created by differentiation with respect to the components of tensor a. 

The GOA consists technically in expanding the logarithm of overlap kernel in powers 
of 7 = a — ck' up to the second order in 7. The coefficients of expansion are the partial 
derivatives of lnX(/3 + ^7, f3 — ^7) at 7 = 0. They are: 



d\nl{(3 + h,(3~h) 



'^ 7-2=---=72=0 

($^(/3)|$(/3))-(<l>(/3)|<l>^(/3)) 



_ (91nX(/3,/3) 

(5.3a) 



{5.3b) 



2($(/3)|$(/3)) 
in virtue of (I5.1a|) and (15.16p . and in accordance with (13. ip . 

_ d'I{f3,f3) 
according to (15.1 aP and (I5.3ap . Hence, the metric tensor defined by (13. ip is 

rw . J^IiMl. (5.4) 

Performing the second order derivative and using additionally (I5.1cp we have finally 

^^'^(«) = ($^(«)|$^(«)). (5.5) 

State |$(q:)) is of the form of (12. 6p . Thus, using (I2.13P we replace the derivatives 
with respect to a^ in (15. 5p by the derivatives with respect to the Euler angles u and 
deformations d. The w-dependence of |$(q;)) is inherent in R{u). Differentiating the 
rotation operator with respect to the Euler angles we obtain the formula 

Lu{uj)R{u) = —R{u)Ju ioTU = x,y,z. (5.6) 

Finally, g'^'^{a.) is obtained in the following form (cf (IB.ip ): 

g^-'{cx) = Y,m^)m^)9ai>{d), (5.7) 

a,b 
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where a, 6 = 0, 2, x, y, z. The intrinsic Cartesian components are: 
gki{d) = {(t)u{d)\(Pi{d)) forfc,/ = 0,2, 
guu{.d) = ^(0(rf)|i^|</>((i)) for M = x,y,^ 
5'A:«(c?) = guu'{d) = for fc = 0,2, u,u' = x,y,z, u^ u', 
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(5.8a) 
{5.8b) 
(5.8c) 



where |0fc((i)) are given by (12. 4p and du is defined in (I2.14p . The intrinsic Cartesian 
components of (15.8 d) vanish owing to the assumed symmetries of \4>{d)) discussed in 
section 12.31 The weight appearing in the Bohr Hamiltonian of (I4.17P is, according to 
(lB.3p . equal to 



VM = \J {gm{.d)g22{,d) - gl^{d))g^^{d)gyy{d)g^^{d). (5.9) 

The intrinsic Cartesian components of matrix f inverse to g and its square root matrix 
n which appear in (I4.18P and (I4.19P can be calculated from (15.8 ap and (15.8&P according 
to flRSJ) and flRGJ) . 

In the Gaussian approximation the energy kernel 'H(a,a') is given by (13. Sp . The 
zero order term in the expansion is 

l-i{OL,OL) 



V[OL] 



v{d) 



{m\H\<P{d)). 



(5.10) 



X(q:, a.) 

The second order term is defined by the matrix of the second order partial derivatives 
of quotient "H/X at 7 = as follows: 

92 



/i'^"(/3) 



( n{i3,(3) 



(5.11) 



where the notation of (I5.3ap is used. Due to (I5.1ap . (I5.3fflp and (15. 4p we have: 



d"H{f3,f3) 



d'jud'j, 



IJ,^ IV 



g^'{l3)v{(3). 



(5.12) 



After replacing the derivatives with respect to 7^ by the derivatives with respect to a^ 
h takes the following form 

h^"'{cx) = - \ {<!>'' {(x)\H\<^'' (ex)) - {<l>{(x)\H\<^^"'{(x)) 

- ($'^(a)|<l>'^(a))($(Q:)|^|$(a)). (5.13) 

Differentiation with respect to a^ by means of (I2.13P gives naturally h in the form of 
( jB.ip with the following intrinsic Cartesian components which do not vanish: 

hkiid) = '^[{Md)\H\Md)) - {<P{d)\H\Md))] - {Md)\Md)) {<P{d)\H\<p{d)) 

for A;,/ = 0,2, (5.14a) 

1 



huuid) 



^dl 



{<P{d)\Hj'Md)) - {<j>{d)\Jl\<p{d)){<p{d)\H\<j>{d)) 
for u = x,y, z. 



(5.14&) 
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The remaining intrinsic components vanish for the symmetry reasons (cf section | 

We see that in order to calculate f , g and h, and then the inverse inertial bitensor 
and the collective potential we should calculate the mean values of the squares of 
angular momenta J^, Hamiltonian H and their products Jl^H within the deformation 
dependent intrinsic ground sate |0(c/)), construct the states \(l)k{d)) and \(t)ki{d)) through 
differentiations of |0((i)) with respect to the deformation parameters d and calculate 
their overlaps and matrix elements of H within them. 

5.2. Inverse inertial functions and potential 

The intrinsic Cartesian components of the inverse inertial bitensor (I4.18P do not depend 
on the Euler angles u and are functions of deformation d only. We refer to the 
corresponding intrinsic components Aah{d)., {a,b = 0,2,x,y, z) of matrix A^i^{cx.) as 
the inverse inertial functions. It results from fl4.18p using (lA.Sp . (IB.ip and (1B.6P that 
the inverse inertial vibrational functions are equal to 

A ij\ _ 922id)hoo{d) - 2gQ2{d)g22id)ho2id) + 5-02 (c^) ^22 (c^) t^ ^K \ 

^"^^^ - {9oo{d)g22{d) - gUd)r ' ^'-'"'^ 

A / n i9ooid)g22id) + 5-^3 (c^)) ^02(0?) - 9o2id)g22id)hooid) - gQ29m{.d)h22{d) ^r, -,,,,. 

^02(0) = 7 TT^ TT^ 2 f^\\2 ' (5.156) 

{goo{d)g22{d) ~ g^2id)y 
A M^ _ 9l2id)hoo{d) - 2go2{d)goo{d)ho2{d) + g^2id)h22{d) /p; 1 p; ^ 

^"^^^ - i9ooid)g22id) - g^dW ^'-''"^ 

and the inverse inertial rotational functions read 

A^d) = %^ (5.15^ 

for u = X, y, z. Hence, the moments of inertia by definition (cf [Ij) are 

/ u^ = ^< = {<P{d)\JMd))' .5 ^g. 

^^ ^ Auuid) {ct>{d)\HJl\ct>{d))-{ct>{d)\PMd)){ct>{d)\H\ct>{d)) 

They resemble the Yoccoz moment of inertia [55] ES] and are in form similar to those of 
a rigid body |3T| . 

Being a quadrupole scalar, the collective potential depends obviously on 
deformation d only. Equation fl4.19p for the potential contains matrices n and a. Their 
intrinsic Cartesian components can be calculated from ( ]B.6p . (IB.SP and (I3.29p . It 



is worth noting that the zero-point energy corrections associated with the rotational 
modes come not only from the rotational components (of type a^u, u = x, y, z) of the 
matrices involved in the corresponding formula for the potential and the derivatives of 
these matrices with respect to d but also from the derivatives with respect to u coming 
from the differentiation of (I2.13P with respect to a^. These latter derivatives can be 
performed using (IB.ip and ( lA.Sp . Obviously, the final result for the potential does not 
depend on the Euler angles. 
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6. Conclusions 

Using tlie GCM we liave generated tlie quadrupole collective excitations from the 
deformation- dependent intrinsic ground state which possesses the D^^ symmetry. To 
consider not only the quadrupole vibrations but also the quadrupole rotations the 
intrinsic state is rotated and the three Euler angles of the rotation are attached to 
the two deformation parameters and form together the real electric quadrupole tensor 
which plays a role of a set of generator coordinates. It turns out convenient to use just 
the quadrupole tensor, the real components of which are the Cartesian coordinates in the 
five- dimensional space, instead of the Euler angles and the deformation parameters. In 
this way, periodic coordinates and a complicated topology of the space are avoided. The 
local GOA has been applied to the overlap and energy kernels. In this case the integral 
Hill- Wheeler equation can be reduced to the differential equation having the form of 
the eigenvalue equation for the Bohr Hamiltonian. The reduction has been performed 
using the Fourier analysis of the energy kernel. The square root of determinant of 
the matrix defining the overlap in the GOA (the metric tensor) is the weight in the 
Bohr Hamiltonian. The exact potential and the inverse inertial bitensor contains the 
exponential function of the five-dimensional Laplacian operator and therefore, can be 
calculated only in an approximation. The simplest approximation consists in replacing 
the exponential operator by the first term of its expansion — the unit operator. This 
approximation corresponds to the potential and the inverse inertial bitensor obtained 
by the usual collective approach to the GCM. The next terms of the expansion can be, 
of course, taken into account in our approach, however, without settling the issue of 
convergence of the procedure. 

When deriving the Bohr Hamiltonian, a transformation of the quadrupole variables 
must be performed, which transforms the metric tensor into the unit matrix. The 
resulting new set of variables is defined well when the square root matrix of the metric 
tensor is an irrotational field (cf jl^)- Obviously, this condition is equivalent to the 
condition which means that the Riemannian space with the metric tensor in question 
being the Euclidean space. The assumption about the flatness of the space is apparently 
made in other papers on the multi-dimensional GOA even if this is not mentioned there 
(cf [T8l [To]). The final form of the Bohr Hamiltonian does not seem to depend on this 
assumption. From the brief report by Kanilah [16] it seems that the assumption is not 
necessary, however, the proof of this fact has not been given there. Is it thus a technical 
condition coming from the fact that we are not able to calculate multi-dimensional 
Gaussian integrals? Or is it an essential condition? A similar problem appears in a 
sense in the case of the Podolsky-Pauli quantization prescription. The form of the 
quantum operator is simply assumed. However, if one would want to derive this form 
from the Schrodinger operator in the Euclidean space, one should assume that the new 
variables are curvilinear variables also in the Euclidean space and the kinetic energy is 
simply proportional to the Laplacian in curvilinear coordinates. It would then mean 
that the classical inertial matrix represents the metric tensor of an Euclidean space. We 
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are left with this open problem. 
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Appendix A. Semi-Cartesian Wigner functions 

The quadrupole tensors transform themselves under rotations by means of the rotation 
matrices or the Wigner functions 'D'i^[u}) which depend on a set of the Euler angles 
defining the rotation [IQ]. It is convenient to introduce the following linear combinations 
of the Wigner functions |41) : 



^moM=^'oM 



1 

V2' 



D,^{u) = ^{V%{u) + Vl_,{uj)), 



D,.{u:) = ^(^'i(c^) + ^'-i(^)), (A.l) 



1 

V2' 



D,y{u) = ^iVl{u)-Vl_,{u)), 



Let us call them the "semi-Cartesian" Wigner functions. The complex conjugate 
functions are 

D^,{u) = D;,{uj) = i-irD^,,{uj). (A.2) 

The orthogonality conditions for the semi- Cartesian Wigner functions take the form: 

2 

Y, m^)D,i{uj) = Ski, (A.3) 

fj.=-2 

J2m^)Duk{co) =5^ (A.4) 

k 

The drift angular momentum operators L„(w), u = x, y, z, act on the semi- 
Cartesian Wigner functions as follows |1T] 

L^{u)D^^{uj) = iV3Df,o{uj) + i/^^2(w), 

L^{u)D^^{uj) = iD^y{uj), L^{uj)Df,2{(^) = -iD^,^{u), 

Ly{u)D^o{u) = iV3D^y{u), Ly{u)Df,^{u) = iD^^i^u), 

Ly{uj)D^y{uj) = -iy/3D^o{u) + iD^2(w), (A.5) 

Ly{uj)D^^{uj) = -iD^^{uj), Ly{uj)D^2iuj) = -iD^y{uj), 
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22 



Lz{u)D^z{uj) 







5 

iD 



m 



,w 



L^{uj)D^y{uj) 






It is seen that all the three components of the drift angular nionientuni change the 
Cartesian indeces of the Wigner functions. 



Appendix B. Symmetric matrices as isotropic functions of the quadrupole 
tensor 

Any symmetric quadrupole bitensor (symmetric matrix 5x5) m^j,(Q;) which is an 
isotropic real function of the quadrupole tensor ck can be determined either by six 
scalar functions oi d = [do, ^2) or by the six intrinsic Cartesian components in the 
following way (cf [Ij): 

,{a) = J2 D^kiuj)D,iiu)mki{d), (B.l) 



m 



iiu\ 



kl 



where the Cartesian indeces k, I run over symbols 0, 2, x, y, z. The Euler angles u 
determine the orientation of the intrinsic system. The Cartesian matrix mki{d) is real 
symmetric and has the following structure: 

/ moo((i) mQ2{d) 

■mo2{d) ■m22{d) 









m{d) 



\ 















myy{d) 











\ 



m. 



(B.2) 



.id) J 



^T^ti^T^xx^T^yy^TT'zz- 



The determinant of m reads 

m = det (m) = (^00^^22 — ''^02)"^^^"^ra"^2 
The Cartesian components of matrix m{cx.) are: 

where C^ = L'^(w = 0). 

If matrix m^j, is not singular i.e. m 7^ matrix (m~^)^j, inverse to m^,y can be 
defined as 



(B.3) 
(B.4) 



m 



flK 



[c.){m-T{c.) 



K- 



(B.5) 



The intrinsic Cartesian matrix inverse to m is: 

m-\d) 

I rn22{d)/ml{d) -mo2{d)/ml{d) 

-mo2{d)/ml{d) moo{d)/ml{d) 





\ 















l/m^^{d) 








l/myy{d) 





1/m, 



\ 



(B.6) 
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If m is positive definite i.e. "moo, ""^22; f^xx-, ''^yyi ^^zz and ml = 771^011122 — "^02 ^^^ 
all positive the positive-definite square root matrix (m^/^)^,^ such that 

m^^{a) = (mV^)^(«)(mi/2)«(c.) (B.7) 

can be defined. It is real symmetric and has the structure of (IB.lj) like m. The 



nonvanishing entries of the intrinsic Cartesian matrix (m^/^)^; are related to components 
iTiki as follows: 

/ 1/2N ^00 + m^ 

f i/2x ^22 + rn^ 



^m jo2 



Vmoo + "^22 + 2m„ 
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